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THE ALGEBRAICAL DEVELOPMENT OF THE ELLIPTIC 
PERTURBATIVE FUNCTION USED IN THE THEORIES 
OF PLANETARY MOTION. 

By R. T. A. Innes, F.R.S.S.Af. 

(Read October 18, 1911.) 

\ 

The methods of computing the motions of the planets devised by 
Lagrange and Laplace and used by the great planetary investigators of 
last century such as Le Yerrier, Hansen, and Newcomb, require an 
expansion of the reciprocal of the distance between two planets. The 
most general expansion is an algebraical one, which, if it is to be useful in 
practical work, requires that the eccentricities and mutual inclinations of 
adjacent planetary orbits should be small—these conditions are fulfilled 
for all the major planets and for the majority of the minor planets of the 
solar system. Hence, great attention has been devoted in the past to this 
algebraical expansion. The most useful, but at the same time the most 
complicated development, is in terms of the mean anomalies of the planets 
concerned, and it is this development which is dealt with here. References 
may be made to the following works :— 

Le Verrier, Les Amiales de VObs. de Paris , vols. i. and x. 

Tisserand, Mecanique Celeste , vols. i. and iv. 

Hill, Collected Works , vol. ii 

Newcomb, American Journal of Mathematics, vol. iii.,and Astronomical 
Papers of the American Ephemeris , vol. v. 

Cowell, Monthly Notices of the Royal Astronomical Society , vol. lxix. 

The object of the present paper is purely practical; it is to put in the 
hands of the computer a development which extends to any order for 
primary or secondary terms, and is substantially ready for use. A 
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numerical application to the theory of the planets Jupiter and Saturn has 
been commenced, but its completion is still far off. 

The most extended tables of the algebraical development of the 
perturbative function in terms of the mean anomalies, are due to 
Le Yerrier (. Les Annales, i. and x.), but a more methodical and general 
development has been given by Newcomb [Astronomical Papers of the 
American Ephemeris, v.). 

Le Verrier’s development has been extended by Boquet to terms of 
the 8th order of both the eccentricities and mutual inclination, and this 
limit has also been adopted by Newcomb. The remarks which follow 
will apply more especially to Newcomb’s development, but with suitable 
changes they can, if necessary, be adapted to Le Verrier’s; it is not, 
however, likely that any one who is familiar with both developments will 

prefer Le Verrier’s—the use of the operator D = aJ^by Newcomb is 

alone sufficient to turn the scale. 

Every term of the development has the form*— 


e n e n’-pyn'' CQS (V ,+jg+j’g’) 


in which k indicates the presence of the factor a 2k (<7 = sine of semi¬ 
inclination). The order of this term will be equal to + 2k. 

ft 

This paper will deal with the properties of the quantities Pthe 
expressions for which form the only difficulty in the development. As it 

o 123 

is very easy to pass from P (henceforth written P) to P, P, P, &c., 
attention will be confined to the former. 

In any given term we find a series of this form— 


eWi?%A-^^,+^A 2 +^f: 2 ^ 2 + &c.) 

(In which the suffices only may be negative and then only under certain 

conditions.) 

The first quantity Pjj! will be called the primary term, those multiplied 
by e 2 or e the secondary terms, those by e 2 e 2 lt e 4 or e\ tertiary terms and so 
on. Terms in which the suffices j or f are negative are secondary or 
higher order terms, although they are closely related to the primary 
terms. 

* It is assumed that the reader is conversant with the notation used by Newcomb in 
the papers above referred to. 


Algebraical Development of the Elliptic Perturbative Function . 303 

But before going on, a further important simplification can be intro¬ 
duced. With Newcomb, we write— 


P^n^A;, 

(Ast. Papers , v., p. 27, &c.) 

in which the operator— 

nr;=n^xn^:. 


Further, as also discovered by Newcomb, the operator II£"! can be 
found from II?- very simply. Hence the real crux is to find expressions 
for llj*. The primary operator in any term will be II", the secondary 
II"+ 2 , the tertiary II" +4 , &c., &c. The expresssions to U] in terms of i and 
D are given by Newcomb. The operators with negative suffices are, as 
already mentioned, at least secondary. 

On account of its complexity, it is almost impossible to give a general 
and useful formula for II" in terms of the D operators. Such a formula 
starts thus— 


(-Iff(n.j)Yl n = + 

+(%*-=<!£V 


+ 1 


| - 2« 2 (« -f) 


D" 

D n ~ 1 


-^(7 - 13?i + 2n 2 + 4/ 2 ) 


In which— 


+ ^_(3?i 3 - 22?i 2 + 45n - 26) 

Jrj 2 (n 2 - 671+4) ++ j. D" - 

&c., &c. 


f(n.j.) — 2"n\ 


n-J 

2 


n.n - l.n - 2 ... n+j + 2 
2 “ 


and f(n.O) = 2 n n\ 


The extent given will merely serve to 7i = 2, although it will always serve 
as a check on the coefficients of the three highest D operators in any value 
of II". 

The primary terms are the more important as the succeeding terms 
are factored by powers of e 2 , &c. In the “ Variation of Elements ” 
method of computing perturbations, the secondary terms become the 
leading terms in several of the coefficients, but in the final expressions for 
the perturbations^of the longitude or mean anomaly the primary terms 
reassert their predominance. Thus in the theory of the motion of the 
earth we find the following perturbations of the longitude 
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Argument. 

Action of— 

Primary. 

Secondary. 

Tertiary. 

I'-X 

Mercury 

+ 0130 

- 0-003 

Nil 

r-\ 

Venus 

+ 3-813 

- 0-002 

Nil 

13 1' - 8X 

Venus 

-1-253 

- 0-003 

Nil 

13 1' - 8X 

Venus 

+ 2-091 

-0012 

Nil 

V-\ 

Mars 

+ 0-102 

+ 0-006 

Nil 

IbV - 8X 

Mars 

+0-161 

- 0 009 

Nil 

21' - 2X 

Jupiter 

+ 1-901 

-0-012 

Nil 


(See for further examples Le Verrier’s “Theory of the Sun ” (Les Annales , 
iv., pp. (9)-(18)). The secondary terms given appear to be the largest 


in this theory.) 

In the action of Saturn on Jupiter we find— 


Argument. 

Primary. 

Secondary. 

Tertiary. 

3l'-2\-w' . 1 

+ 26-79 

- 010 

+ 002 

Largest individual term of 

the Great Inequality . 

-1522-78 

+ 11-18 

+ 0-41 


(See Les Annales , xi., p. 121, &c.) 

In the last case the next term may be estimated at about 0”-04. For 
such a quantity to have any real significance we must know the mass of 
Saturn correctly within its - 30^00 part, and the annual motion of Jupiter 
within 0"*00008. Both of these requirements are beyond our power of 
attainment; in short, terms of the tertiary order are of very little import¬ 
ance, whilst those of the next order are quite negligible in the theories to 
which an algebraical development can be applied. 

It may, however, be necessary to include such terms in the considera¬ 
tion of the secular perturbations, but in this case it is possible to give 
general expressions of great simplicity, extending to any powers of the 
eccentricities and mutual inclinations, and these will be found in the last 
table of the present paper. 

Chessin has shown {Ast. Journal , No. 442, 1898) that it is possible to 
simplify the expressions for the computation of the values of the Newcomb 
operators ; he gives the following example :— 

Newcomb— 

768011$ = D 6 + {Si - 25)D5 + (20* 2 -130*' + 185)D4 + (- 80^ 2 + 340^ - 255)Ds 
+ (_ 80i4 + 640*3 -1765* 2 + 2280* - 1466)D 2 
+ (-128*5 +1360*4 - 5280*3 + 9535* 2 - 8454* + 3096)D 
+ (- 64* 6 + 800*5 - 3740*4 + 8200*3 - 8588* 2 + 3608*). 






















Algebraical Development of the Elliptic Perturbative Function. 305 
’Chessin— 

n*=6n* + 2 (i - 5)n | - (f i - 2)nj- (i - i)n 3 - *(2i- i)nj - T h(5i - i)n; - T v‘ 

Chessin’s formula has two advantages; firstly, it is much shorter p 
secondly, it throws the computation of the secondary terms on to those of 
the first. In ordinary course the primary terms must be computed first, 
and it is fitting that they should form the basis for the secondary and 
higher terms. But neither Newcomb nor Chessin gave a general 
expression for the operators, and, in point of fact, it requires a very 
lengthy calculation to check, far less to extend, Newcomb’s tables. 

Let us suppose that we are dealing with the primary operator n*; 
Newcomb’s formula starts thus— 

10321920n^ = 256i 8 - 8960^+130592i 6 -1023120*5+ 462807** + &c. 
plus terms in D, D 2 ... to D 8 . 

The figures bear no obvious relation to each other, and have to be 
taken on trust; besides this, the computation for a particular value of i is 
exceedingly arduous. The length of the expression for n^, n;°, &c., is 
appalling ; it is true that the use of the Chessin formulae would shorten 
the calculation, but Chessin does not give the constants his method 
requires for operators beyond and n*. Chessin’s constants for II” are, 
however, not difficult to find, as they are connected with certain Besselian 
series for e used in Hansen’s planetary theory. 

Newcomb gives, on p. 13 of vol. v. of the Astronom . Papers of the 
American Ephemeris , the general recurrent formula by which he builds 
up II”; it is— 

2(n + l)n”£ = (Kp + fe;D)n; 

+ (K / i + KD) u n n z\ 

+ . 

+ (kRlfi + h#\D) n° 

in which the factors k and h depend upon the developments of the 
equation of the centre and the logarithm of the radius vector in terms 
of the mean anomaly. The succession of the values of the k’ s given by 
Newcomb (p. 22 loc. cit .) are— 

9 5 13 103 1097 1223 47273 556403 
2’ 4’ 24’ 192"’ 160' 4608* ~40320 

What the next fraction will be is by no means clear, but it will be seen 
immediately that the values are really not wanted, as a transformation at 
once discloses the values of the successive terms. Let us write out the 
Newcomb values for n;, II^ in the following form :— 
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11°. 

2n;. 

8n). 

4811). 

-(2t + D) 

-l 


( 

— (hi + 3D) 

-(2* + D) 

-1 

( 

-(26+ 17D) 

— (10* + 6D) 

-(2* + D) 

-1 i 

&c. 

&e. 




(a) 

(*) 

(e) 


which is read as follows :— 

0= -(2* + D) n°-2n; &c. 

Subtract now 3(a) from (b), and generally subtract whole multiples of 
(a) ( b )... from the following letter and the following table results:— 


11°. 

2n;. 

8n;. 

4811). 

- (2i-f D) 

i 

-l 

- 2f + 3 - D 

-1 


2 i 

2f -1 

- 2i -f 6 - D 

-1 1 

&c. 

&c. 




(ay 

(by 

(cy 


The final results are presented in Table I.a } which contains precepts for 
its indefinite extension. The n°;$ will be found in Table Lb. It will also- 
be seen that the property— 

falls naturally into the scheme of these tables. 

In practice we first compute— 


and— 

then- 


2n;, 8n 2 2 , 48n 3 ... &c. 
2n° : ;, 8n° ; 2 , 48ns ••• &c. 
4n; : ;, i6n 2 ;;, 96n 3 ;; ... &c. 


and so on, step by step, each individual step being simple enough, and 
natural checks on the results abound. This method requires, however, 
the quantities D 5 n° where q- 1 is equal to the highest power of e or e iy 
which has to be included. Thus in actual work we must compute— 

2n;, 8n 2 , 48n 3 ... &c. 

D2n;, D8n 2 2 , D48n 3 ... &c. 

D 2 2n;, D 2 8n 2 , &c. ... &c. 
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The presence of the D operator in no way alters the tables, thus— 

0 - (-2i-D)ns-2n:. 

0 = D*( - 2 i - D)n° - D ? 2n;. 

And so on. 

As pointed out in the tables, the forms for or are included. 

It must be remembered that such an operator as n_[ cannot act alone ; in 
actual work it will be joined to another of at least the same order but 
with a positive suffice, the final form being R_\\\ or and so on. In fact, 
although of primary form, such are closely connected with the secondary 
terms. As the secondary, tertiary, and higher terms are seldom required, 
it is, as already mentioned, advisable to make them dependent on the first 
order terms. In the tertiary terms we first meet with the quantity i 2 , but 
it can be eliminated by the introduction of an 7D operator which will 
make the numerical work simpler. 

The tables for the secondary terms (Tables II.A and II.b and III.A 
and III.B) require no further explanation—the examples attached to each 
show their use. 

The paper concludes with a symbolic expression for the non-periodic 
or secular part of the perturbative function, which by simple multiplica¬ 
tions will give the expansion to any powers of the excentricities and 
mutual inclination. Examples are added, which can be compared with 
the known results given by Newcomb and other investigators. 

The happy introduction by Newcomb of the two operators—one of 
which is purely symbolical—D and IT, has vastly simplified the develop¬ 
ment of the perturbative function. The curious properties of the D 
operator have been dealt with in the Monthly Notices of the Eoyal 
Astronomical Society for June and December, 1909. 


Johannesburg, 

September 5, 1911. 
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This example is chosen because it shows how many of the higher terms can be independently calculated in 
duplicate. 
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Examples. 96n_|;;’= — 204811°;' - QGIIIJ. (See for this example, also Table II.A.) 
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Non-Periodic Part of the Elliptic Perturbative Function. 
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Third line. To each angle in the first line both add and subtract 2^ — 2^!, &c., &c. 
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